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We calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova sim- 
ulations spanning a range of massive stars. Most of the simulations are carried out to times late 
enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode 
of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inex- 
orably rises as the proto-neutron star (PNS) core shrinks. We demonstrate that the GW emission 
is excited mostly by accretion plumes onto the PNS that energize modal oscillations and also high- 
frequency (“haze”) emission correlated with the phase of violent accretion. The duration of the 
major phase of emission varies with exploding progenitor and there is a strong correlation between 
the total GW energy radiated and the compactness of the progenitor. Moreover, the total GW 
emissions vary by as much as three orders of magnitude from star to star. For black-hole formation, 
the GW signal tapers off slowly and does not manifest the haze seen for the exploding models. For 
such failed models, we also witness the emergence of a spiral shock motion that modulates the GW 
emission at a frequency near ~100 Hertz that slowly increases as the stalled shock sinks. We find 
significant angular anisotropy of both the high- and low-frequency (memory) GW emissions, though 


the latter have very little power. 


The theory of core-collapse supernova (CCSN) explo- 
sions has been developed over the last six decades and 
is now a mature field at the interface of gravitational, 
particle, nuclear, statistical, and numerical physics. The 
majority of explosions are thought to be driven by neu- 
trino heating behind a shock wave formed upon the colli- 
sion of the rebounding inner core with the infalling man- 
tle of the Chandrasekhar mass birthed in the center of 
stars more massive than ~8 Mo[1-5]. After implosion 
ensues, this inner white dwarf core, with a mass near 
~ 1.5Mo and a radius of only a few thousand kilome- 
ters, requires only hundreds of milliseconds of implosion 
to achieve a central density above that of the atomic nu- 
cleus. At this point, the inner core stiffens, rebounds, and 
collides with the outer core, thereby generating a shock 
wave that should be the supernova explosion in its in- 
fancy. However, detailed 3D simulations [6-22] and phys- 
ical understanding dictate that this shock generally stalls 
into accretion, but is often reenergized into explosion by 
heating via the neutrinos emerging from the hot, dense, 
accreting proto-neutron star (PNS), aided by the effects 
of vigorous neutrino-driven turbulent convection behind 
that shock [15, 22-26]. The delay to explosion can last 
another few hundred milliseconds, after which the explo- 
sion is driven to an asymptotic state in a period of from 
~a few to ~10 seconds. An extended period of neutrino 
heating seems required [5, 15, 27]. The shock wave then 
takes a minute to a day to emerge from the massive star, 
and this emergence inaugurates the brilliant electromag- 
netic display that is the supernova. The outcomes and 
timescales depend upon the progenitor core density and 
thermal structure at the time of collapse [28, 29], which 
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itself is an important function of progenitor mass, metal- 
licity, and rotational profile. If a black hole eventually 
forms, the core must still go through the PNS stage, and 
it is still possible to launch an explosion, even when a 
black hole is the residue. There is never a direct collapse 
to a black hole. A small fraction of supernova (hyper- 
novae?) may be driven by magnetic jets from the PNS 
if the cores are rotating at millisecond periods. Other- 
wise, magnetic effects are generally subdominant, but of 
persistent interest in the context of pulsar and magnetar 
birth [30-36]. 

Though this scenario is buttressed by extensive simu- 
lation and theory, and most 3D (and 2D) models now ex- 
plode without artifice [8, 12, 13, 15, 17, 20-22, 26, 37, 38], 
direct verification of the details and the timeline articu- 
lated above are difficult to come by. However, the neu- 
trino and gravitational-wave signatures of this dynamical 
event would allow one to follow the theoretically expected 
sequence of events in real time. The detection of 19 neu- 
trinos from SN1987A was a landmark [39, 40], but little 
was learned, other than that a copious burst of neutrinos, 
whose properties are roughly in line with theory, attends 
core-collapse supernovae and the birth of a compact ob- 
ject. The real-time witnessing of the events described as 
they unfold is the promise of gravitational-wave (GW) 
detection from a supernova explosion. The detection of 
these GWs and the simultaneous detection of the neutri- 
nos overlapping in time is the holy grail of the discipline. 


The CCSN gravitational wave signal [41-54] from 
bounce through supernova explosion and into the late- 
time (many-second) proto-neutron star (PNS) [32, 55, 56] 
cooling phase (or black hole formation) can be de- 
composed into various stages with characteristic fea- 
tures, frequency spectra, strains and polarizations [57]. 
GWs are generated immediately around bounce in core- 


collapse supernovae by time-dependent rotational flat- 
tening (if initially rotating, see e.g. [58]) and prompt 
post-shock overturning convection (lasting tens of mil- 
liseconds) [24, 59], then during a low phase (lasting 
perhaps ~50—200 milliseconds) during which post-shock 
neutrino-driven turbulence builds, followed by vigorous 
accretion-plume-energized PNS modal oscillations (pre- 
dominantly a mixed f/g-mode early, then a pure {mode 
later, see §I1A) [41, 51]. If an explosion ensues, these 
components are accompanied by low-frequency (~1—25 
Hertz (Hz)) gravitational-wave “memory” (see §II F) due 
to asymmetric emission of neutrinos and aspherical ex- 
plosive mass motions [60-63]. The duration of the more 
quiescent phase between prompt overturning convection 
and vigorous turbulence depends upon the seed pertur- 
bations in the progenitor core [15], of which the duration 
of this phase is diagnostic. Relevant for the GW sig- 
nature is the equation of state (EOS) [59, 64-66], the 
rate of neutronization and neutrino cooling of the core 
(56, 67], and the stellar core’s initial angular momentum 
and mass density distributions. Measurable GW signa- 
tures of rotation require particular progenitors that ro- 
tate fast, while all other phases/phenomena are expected 
to be operative in any core-collapse context, except mat- 
ter memory, which requires an asymmetric explosion. 
The rotational signature is primarily through its depen- 
dence upon the ratio of rotational to gravitational energy 
(T/W) [18, 44, 50], and, for fast rotating cores, to the de- 
gree of initial differential rotation. Black hole formation 
will approximately recapitulate the sequence followed by 
neutron star formation, except that when the black hole 
forms long after collapse the GW signal ceases abruptly 
[68] and that due to the shrinking stalled shock radius a 
spiral mode is excited that quasi-periodically modulates 
its GW signal (see §II B). The termination of the neutrino 
emission is correspondingly abrupt. Though these phases 
are generic, their duration, strain magnitudes, and degree 
of stochastic variation and episodic bursting (due in part 
to episodic accretion and fallback) is a function of mas- 
sive star progenitor density structure, degree of rotation, 
and the chaoticity of the turbulence. 


Though neutrino-driven convection generally over- 
whelms manifestations of the “SASI” (Standing Accre- 
tion Shock Instability [69, 70]), the SASI is sometimes 
discernible as a near 100—200 Hz subdominant compo- 
nent. However, when the explosion aborts and the av- 
erage shock radius sinks deeper below ~100 km, the so- 
called “spiral SASI” [46, 71, 72] emerges (§II B). This spi- 
ral rotational mode has a frequency of ~100-200 Hz, is 
interestingly polarized, and if present clearly modulates 
both the neutrino and GW signatures until the general- 
relativistic instability that leads to black hole formation. 


Therefore, each phase of a supernova has a range of 
characteristic signatures in gravitational waves (GWs) 
that can provide diagnostic constraints on the evolution 
and physical parameters of a CCSN and on the dynamics 
of the nascent PNS. Core bounce and rotation, the exci- 
tation of core oscillatory modes, neutrino-driven convec- 


tion, explosion onset, explosion asymmetries, the mag- 
nitude and geometry of mass accretion, and black hole 
formation all have unique signatures that, if measured, 
would speak volumes about the supernova phenomenon 
in real time. 


For this paper, we have run a broad set of 3D simu- 
lations (11 progenitors in total, with progenitors ranging 
from 9 to 23 Ma) out to the late post-bounce times (up 
to ~6 seconds post-bounce). These late-time 3D simula- 
tions illustrate the sustained multi-second gravitational 
wave signal and are the longest 3D core-collapse simula- 
tions with sophisticated neutrino transport performed to 
date. Simulations out to one second or less do not capture 
the entire time evolution of the gravitational wave signal. 
We find that the gravitational wave signal persists out to 
late times for all models, is strongly correlated with the 
turbulence interior to the shock (see §II G), and is not 
correlated with proto-neutron star convection (§IIG 1). 
We see a memory signature at <25 Hz associated with 
large-scale ejecta for our exploding models, and a spiral 
SASI signature at ~100 Hz for the non-exploding models. 
All our models show an early prompt-convection phase at 
~50 milliseconds (ms) associated with a negative entropy 
gradient interior to the stalled shock front at ~150 km. 
For most models, this is followed by a quiescent phase of 
duration <50 ms, after which the strain grows in associ- 
ation with turbulent motion interior to the stalled shock. 
The 23-M© model shows an interesting exception, as the 
strain illustrates another feature lasting from ~175 to 
350 ms, coincident with the shock receding before reviv- 
ing. As the accretion rate and turbulence diminish, so 
too does the strain. However, at late times, we see a 
consistent offset in the strain associated with matter and 
neutrino memory in the exploding models (§I F). 


We now proceed to a more detailed discussion of our 
new results that span a broad progenitor mass range, 
capture for the first time and for most models the entire 
gravitational wave signal of CCSNe, and do not suffer 
from Nyquist sampling problems. In this paper, we fo- 
cus on initially non-rotating progenitors, whose general 
behavior should also encompass models for slowly rotat- 
ing initial cores. We note that even initially non-rotating 
3D CCSN models experience core spin up due to stochas- 
tic fallback [73], and this effect is a normal by product 
of sophisticated 3D simulations. In §I, we provide in- 
formation on the simulation suite and its characteristics 
and briefly describe the various models’ hydrodynamic 
developments. Then, in §II we present our comprehen- 
sive set of findings concerning the complete gravitational- 
wave signature of initially non-rotating core-collapse su- 
pernovae, partitioned into subsections that each focus on 
a different aspect of this signature and its import. In 
subsection gII A, we lay out the basic signal behaviors as 
a function of progenitor. This section contains our ma- 
jor results and then is followed in §IIB by a digression 
into the GW signature of black hole formation. In §IIC, 
we present an interesting finding concerning the depen- 
dence of the total radiated GW energy on compactness 


[74] and in §II D we note the avoided crossing that is uni- 
versally manifest in all CCSN GW spectra and seems a 
consequence of the presence of inner proto-neutron star 
(PNS) convection and its growth with time. In §IIE we 
discuss the solid-angle dependence of the matter-sourced 
GW emission, both at high and low (matter “memory” ) 
frequencies, and in §II F we present our results concern- 
ing neutrino memory at low frequencies. Then, we tran- 
sition in §IIG to a discussion of the predominant exci- 
tation mechanism. Finally, in §III we recapitulate our 
basic findings and wrap up with some observations. 


I. SETUP AND HYDRODYNAMICS SUMMARY 


We present in this paper a theoretical study of the 
gravitational-wave emission of eleven core-collapse super- 
novae (CCSNe) progenitors, from 9 to 23 Mo, evolved in 
three dimensions using the radiation-hydrodynamic code 
FORNAX [75]. To calculate the quadrupole tensor and 
the gravitational-wave strains we employ the formalisms 
of references [76] and [77] (see also §A) and we dump 
these data at high cadences near the LIGO sampling rate 
(Table I). The progenitor models were selected from [78] 
for the 14-, 15.01-, and 23-Mo models and from [79] for 
models between 9 and 12.25 Mọ. The radial extent of 
the models spans 20,000 kilometers (km) to 100,000 km, 
generally increasing with progenitor mass. All the mod- 
els (except the 9-M model on Blue Waters, which had 
648 radial zones) were run with 1024x 128x256 cells in 
radius, 0, and ¢. We employ 12 neutrino energy groups 
for each of the ve, De, and “v,,”s followed (see [21, 73]) 
and the SFHo equation of state [80]. The progenitor 
models are non-rotating, though some degree of rotation 
is naturally induced due to fallback [73]. These simula- 
tions include two of the longest 3D CCSNe simulations 
run to date, a 11-Mọ model evolved past 4.5 seconds 
post-bounce, and a 23-Mọ model evolved to ~6.2 sec- 
onds post-bounce. All of our models explode except the 
12.25- and 14-Mọ progenitors, and all models besides 
the 14-Mo are evolved beyond one second post-bounce. 
We include four simulations of the 9-Mọ model on dif- 
ferent high-performance clusters (Frontera, Theta, and 
Blue Waters) at various stages of code evolution. The 
low-mass 9-Mọ models asymptote early on in diagnostic 
quantities such as explosion energy and residual mass, 
and one iteration has been evolved past two seconds post- 
bounce. Our models, run times, and explosion outcome 
are summarized in TableI. Several of the models have 
been published before, including three of the four 9-Mọo 
iterations (the fourth, the longest simulation, is new) in 
references [81], [21], and [73]. 

In Figurel, we plot the density profiles against en- 
closed mass for all eleven models studied here. Note the 
association of the silicon-oxygen (Si/O) interface density 
drop (see, e.g. [5, 14, 21, 28, 29, 38, 82, 83]) with the 
onset of successful shock revival. The 14-Mọ and 12.25- 
Mo progenitors lack such a strong interface and do not 


explode. Low-mass progenitors (e.g. 9-9.5 Mo) have a 
steep density profile and explode easily (though still with 
the aid of turbulent convection). Models 11-, 15.01-, and 
23-Mo have Si/O interfaces successively further out, and 
explode successively later. 

We show the angle-averaged shock radii at early and 
late-times in Figure 2. All models except the 14 and 12.25 
Mo models explode, with an approximate correlation be- 
tween progenitor compactness [74] and explosion time. 
The two non-exploding models experience ~10 ms oscil- 
lations in the shock radii due to a spiral SASI that man- 
ifests itself after ~350 ms. We also see a longer secular 
timescale oscillation of ~70 ms in the 12.25-Mọ and 14- 
Mg black-hole formers. We summarize the eleven models 
in Table I. The more massive 15.01- and 23-Mo progeni- 
tors explode later, with the latter showing shock revival 
only after ~0.5 seconds post-bounce. A later shock re- 
vival time, again ~0.5 s, was also seen for the 25-Mo 
progenitor in [21]. After the first ~500 ms, the shock ve- 
locities settle into approximately asymptotic values that 
range from 7000 to 16000 km s~!, inversely correlated 
crudely with the progenitor mass (see also [84]). 


II. RESULTS 


A. General Gravitational-Wave Signal Systematics 
of Core-Collapse Supernovae 


As stated, we highlight for this study of the GW sig- 
natures of core-collapse supernovae eleven of our recent 
initially non-rotating 3D FORNAX simulations. Care has 
been taken to calculate the quadrupole tensor with a high 
enough cadence to avoid Nyquist sampling problems, and 
we have been able to simulate to late enough times to cap- 
ture what is effectively the entire GW signal after bounce 
for a large subset of the models. Table I provides the du- 
ration of each simulation and the minimum Nyquist fre- 
quency achieved during each run. We will discuss both 
matter and neutrino contributions to gravitational wave 
energy, and the relevant equations are summarized in Ap- 
pendix A. 

In Figure 3, we plot the plus (black) and cross (blue) 
polarizations in the x-direction of our computational grid 
of the strain multiplied by the distance to the source. 
Other orientations yield qualitatively similar numbers for 
the higher frequency components that dominate the GW 
power. However, there is a large variation at low frequen- 
cies (<25 Hz) of the matter and neutrino memories with 
solid angle (see §II F). Note that the x- and y-axes cover 
different ranges for each model. A red star on the pan- 
els indicates the rough time of explosion, defined loosely. 
As Figure 3 demonstrates, all the exploding models tran- 
sition through similar phases. During the first ~50 ms 
there is a burst of emission due to prompt overturning 
convection driven by the negative entropy gradient pro- 
duced behind the shock wave as it stalls. The detailed 
time behavior of this overturn will depend on the initial 


accreted perturbations, which will set the number of e- 
folds to the non-linear phase. However, the basic behav- 
ior and timescales are broadly similar. Figure 4 focuses 
on this early first 0.25 seconds. For the lower-mass pro- 
genitors, explosion (the red star) ensues towards the end, 
or not long after, the prompt signal, and this is followed 
by the early growth of the second phase. For the more 
massive progenitors (such as the 23-Mọ model), the on- 
set of explosion can be much later. As Figure 3 shows, 
the growth phase of the GW emission continues beyond 
what is shown in Figure 4 to a strong peak. That peak 
phase is powered by the accretion of the infalling plumes 
during explosion. A core aspect of 3D core-collapse ex- 
plosions is the breaking of spherical symmetry that al- 
lows simultaneous accretion in one direction and explo- 
sion in another [5, 21, 28, 72]. For exploding models, the 
infalling plumes that strike the surface of the PNS can 
achieve supersonic speeds before impact. For the black 
hole formers (the 14-Mo and 12.25 Me models here), the 
accretion is maintained, but impinges upon the PNS core 
subsonically. This will have interesting consequences we 
discuss in §II B. 

Figures 3 show that the lower-mass models have 
smaller strains and that the phase of high strain lasts for 
a shorter time. For the 9-Mo through 9.5-Mọ models, 
much of the GW emission subsides by ~0.25-0.5 seconds, 
while the high phase lasts ~1.2 seconds for the 23-Mo 
model and continues beyond ~1.0 and ~1.5 seconds for 
the 11-Mo and 15.01-Mọ models, respectively. These dif- 
ferences reflect the differences in the initial density pro- 
files (Figure 1) and the compactness (see also Figure 11). 

After this vigorous phase, the pounding of the accre- 
tion plumes subsides, but the signal continues at a low 
amplitude. Though as much as ~95% of the GW energy 
emission has already occurred, the f-mode continues to 
the latest times we have simulated as a low hum of pro- 
gressively increasing frequency". Hence, we see univer- 
sally for the exploding models a transition from a high- 
amplitude, lower-frequency stage (<0.3-1.5 seconds, de- 
pending upon the progenitor) to a lower-amplitude high- 
frequency stage (>1.5 seconds). As Figure 3 indicates, 
for the exploding models a very-low frequency memory 
is superposed that represents a permanent metric strain. 
There is no such matter memory signal for the black- 
hole formers (§IIF), but the accretion phase continues 
for them to very late times, abating only slowly as the 
mantle continues to accrete the mass of the outer man- 
tle until the general-relativistic instability that leads to 
a black hole ensues. 

It has too often been thought that strain signals such 
as are depicted in Figures 3 and 4 are too noisy to be tem- 
plated cleanly, and this to a degree is true. There is a lot 
of stochasticity due to chaotic turbulence. However, the 
frequency content of these signals tells a different story. 


1 Sonifications of the signals are available upon request. 


In Figures5 and 6, we plot spectrograms of the gravi- 
tational wave power versus time after bounce for our 3D 
models and see distinct structures. The most obvious fea- 
ture is the f/g-mode [41, 51, 85-87] from ~400 Hz early, 
rising to ~1000-3000 Hz after ~0.8 seconds after bounce. 
It is in this band that most of the emitted GW power of 
supernovae resides (see Table I). This is a natural conse- 
quence of the fact that the peak in the eigenfunction of 
the f-mode is in the PNS periphery where the collisions of 
the accreta with the core are occurring. Hence, the exci- 
tation and the fundamental f-mode eigenfunction overlap 
nearly optimally. Associated with this feature in the ear- 
lier phases is a dark band near ~1000-1300 Hz. This 
has been interpreted as a manifestation of an avoided 
crossing [51] between a trapped £ = 2 g-mode and the 
€ = 2 fmode. All the spectrograms for all our mod- 
els show the same modal interaction, though at slightly 
different frequencies. For instance, at two seconds after 
bounce, the f-mode frequency is ~1.75 kHz, 1.8 kHz, 2 
kHz, and 2.5 kHz for the 9-, 11-, 12.25-, and 23-Mọ mod- 
els, respectively, reflecting the variation in model PNS 
masses. Early on power is in the lower frequency compo- 
nent (mostly a trapped g-mode, mixed with the fmode), 
and then it jumps to the higher frequency component 
(mostly the f-mode). This modal repulsion, or “bump- 
ing,” is a common feature in asteroseismology [88] and 
seems generic in core-collapse seismology. 


All the models show the early prompt convection 
phase, with power from ~300 to ~2000 Hz. After this, 
all the models manifest a “haze” of emission that ex- 
tends above the f/g-mode to frequencies up to ~2000 to 
~5000 Hz. The duration of this haze is from ~0.25 to 
~1.5 seconds and tracks the phase of vigorous accretion 
(see §II G). Individual PNS pulsation modes, likely the 
£L = 2;n = 1 p-mode, can also be seen superposed in this 
“haze” and extending beyond it to later times. This is 
particularly the case for the 12.25-, 15.01-, and 23-Mo 
models. It is only for the models with the most vigorous 
accretion onto the core that this mode is clearly seen at 
later times. 


The origin of this haze is still a bit unclear, though it is 
definitely excited by the pummeling accreta (§II G). The 
p-mode frequencies of the PNS for radial node numbers 
from 1 to 10 reside in this space and we could be seeing an 
overlapping and unresolved superposition of these modes. 
However, exploding models experience simultaneous ex- 
plosion and accretion, the latter through infalling funnels 
that are few in number, can achieve supersonic speeds, 
and dance over the PNS surface. It is likely that the 
time-changing quadrupole moment of these funnels as 
they impinge upon the PNS surface is the source of this 
power. The timescales of their deceleration are about 
right, those timescales have a spread which could trans- 
late into a broad feature, and at any particular time they 
represent a low angular order perturbation. Importantly, 
however, we don’t see this haze for the black-hole formers 
12.25-M. and 14-Mo. It is only for the exploding mod- 
els that there is a breaking of spherical symmetry that 


results in simultaneous explosion and supersonic funnel 
infall. Figure 7 portrays two snapshots of the Mach- 
number distribution of the inner 100 km of the residue of 
the 23-Mọ model, clearly showing such funnel collisions. 
The accretion onto the proto-black-hole is always sub- 
sonic. Though the haze constitutes at most only a few 
10’s of percent of the total emission, its origin is clearly 
an interesting topic for future scrutiny. 

In Figure8, we plot spectrograms for two represen- 
tative models of the effective strain versus time after 
bounce. The effective strain is defined as the average 
of both strain polarizations, 


heg = 0.5 (h4 + hy) i (1) 


This figure provides a focus on the low-frequency regions. 
We see for the 14-Moọ black-hole former some power near 
~100 Hz, which we identify with the spiral-SASI [71] 
(see §II B). Such a mode emerges in our model set only 
for the proto-black-holes (see §II B). In the red blot in 
the lower left-hand corner of the left panel may be a 
signature of the traditional SASI [69]. We generally see 
little evidence in the GW signature of this SASI, but 
always see the spiral SASI when the explosion is aborted 
and the stalled shock recedes. In the right panel of Figure 
8, the red band is a signature of the matter memory 
associated with the asymmetrical explosion of the 23- 
Mo model. Whether the traditional SASI is seen to the 
left of this band is unclear, but the early recession of its 
shock before explosion might be conducive to its brief 
appearance. 


B. Signatures of Failed Explosions - The Prelude 
to Black Hole Formation 


If a black hole forms by late-time fallback after many, 
many seconds to hours after the launching of a stalled 
shock that seemed to herald a successful explosion (but 
didn’t), then the GW signal will be similar to those seen 
in the context of successful explosions. If, however, the 
stalled shock is never “reignited,” it will slowly settle to 
progressively smaller radii and the mantle of the pro- 
genitor core will continue to accrete through it onto the 
PNS. Eventually, the fattening PNS will experience the 
general-relativistic instability to a black hole, at which 
time the GW emission abruptly terminates within less 
than a millisecond. This latter phase could take many 
seconds to many minutes to reach?. The GW signature 
of this modality of black hole formation, representatives 
of which are our 14- and 12.25-Mo models, has particular 
diagnostic features that set these evolutions apart. First, 
the breaking of symmetry that results in the simultane- 
ous accretion of lower-entropy plumes with the explosion 


2 For the 14-Mọ model, we estimate a black hole formation time 
(using a maximum baryon mass of 2.477 Mọ from Steiner et al. 
[80] at the onset of collapse) of ~500 seconds. 


of high-entropy bubbles does not occur. The result is 
that for this channel of black hole formation the infalling 
plumes do not dance over the PNS, do have high Mach 
numbers, and don’t excite the higher-frequency “haze” 
that we have identified for the exploding models seen in 
the associated spectrograms (Figures 5 and 6). We do 
see in Figure 6 power not only in the dominant f-mode, 
but weakly in an overtone p-mode as well. However, as 
shown in Table I for the 12.25 Mo black-hole former, the 
fraction of the total GW energy radiated in the f-mode 
is correspondingly higher, as much as ~95% of the total, 
than for the exploding models that also generate power 
in the haze. 

This channel of black hole formation also experiences 
the emergence of what we identify as the spiral-SASI 
[71]. This is seen in Figure 2 in the clear ~100-200 
Hz periodicity of the late-time mean shock position of 
both the 14- and 12.25-Mọ models after ~300—400 mil- 
liseconds after bounce and very clearly in the spectro- 
gram of the shock dipole depicted in Figure 9. Gener- 
ally, this feature emerges after the mean stalled shock 
radius sinks below ~100 km and is not seen in exploding 
models. The timescale of the periodicity scales roughly 
with AR,/cs + ARs/Vacc, where Rs, Cs, and Vace are the 
shock radius, speed of sound, and post-shock accretion 
speed [70]. 

Another feature seen most clearly in Figure 2 in the 
context of these black-hole formers is a much longer- 
timescale modulation of the mean shock position with a 
period near ~70 ms. Not seen clearly in the GW spectro- 
grams or strain plots (though there may be a hint in the 
strain plot for the 14-Ms model), this oscillation may 
be due to a global pulsation mode associated with the 
neutrino heating, cooling, and transport of the mantle, 
but this speculation remains to be verified. Nevertheless, 
this feature has never before been identified in studies of 
3D CCSNe and is interesting in itself. Finally, as Fig- 
ure 3 suggests for the 14- and 12.25-Mọ models, since 
those cores that form black holes by this channel do not 
explode, they are expected to have no net low-frequency 
matter memory component. 


C. Total Gravitational Wave Energy Radiated 


In Figure 10, we plot versus time the integrated radi- 
ated gravitational wave energy due to matter motions. 
Model 23-Mo radiates the most gravitational wave en- 
ergy (~3.0x10*6 erg, or ~2x1078 Mo c? after ~5 sec- 
onds), while the collection of 9-Mọ models radiates the 
least. There are a few important features of this plot. 
First, we see that we have captured what is basically the 
entire GW signal for many of the models (the 15.01-Mo, 
14-Mo, and 12.25-Mo model emissions are still climb- 
ing). Within ~1.5 seconds, most exploding models have 
radiated >95% of the total energy to be radiated, and 
after ~2 seconds they have radiated >98%. Table I pro- 
vides the total energy radiated via the f/g-mode, as well 


as the fraction of this total radiated in the f-mode after 
1.5 seconds. Not shown in Table I is the fact that more 
than 95% of the total GW energy radiated after 1.5 sec- 
onds is via the f-mode. As indicated for the 12.25-Mo 
model in Table I and Figure 10, due to continuing accre- 
tion the black hole formers radiate to later times than 
the exploding models, and this mostly in the f-mode. 
Figure 10 shows that the various phases described in 
SI A are recapitulated via stair steps until finally asymp- 
toting. Moreover, the continuum of models highlighted 
in this paper demonstrate collectively that the radiated 
GW signal energies vary by as much as three orders of 
magnitude from the lowest-mass to the higher-mass mod- 
els. This is a consequence of the differences in their initial 
density profiles (see Figure 1), and directly from the re- 
sulting mass accretion histories. Even more directly, as 
Figure 11 demonstrates, there is a strong monotonic rela- 
tion for exploding models between the total GW energy 
radiated and the compactness [74] of the initial progen- 
itor “Chandrasekhar” core”. Though compactness does 
not correlate with explodability [5, 21, 28, 29, 83], it does 
seem to correlate with residual neutron star mass, radi- 
ated neutrino energy, and, as now indicated in Figure 11, 
the total gravitational-wave energy radiated. In fact, we 
derive a power-law with index ~0.73 between the two. 
Finally, we note that the collection of 9-Mọ models 
don’t behave exactly the same. This is due to the fact 
that these models were simulated with slightly different 
code variants (as we continued to update and upgrade 
FORNAX); one (the 9a model) had small artificial im- 
posed perturbations in the initial model and three dif- 
ferent supercomputers were used. The natural chaos in 
the flow and in the simulations will pick up on any slight 
variations and amplify them, with the result that the evo- 
lution can be slightly different. Though the qualitative 
behavior of all four of these 9-Mọ models is the same, 
they exploded at slightly different times (see Figure 3), 
with the resulting different developments of their GW 
signals. As a consequence, the radiated energy varies by 
about a factor of two. In a crude sense, we might view 
this spread as an imperfect indicator of the likely spread 
in Nature due to the chaos of turbulence for the “same” 
progenitor, but we have certainly not demonstrated this. 


D. Avoided Crossing and Trapped g-mode 


A distinctive feature seen clearly in the spectrograms 
for all the models (see Figures 5, 6) is a dark band near 
~1000—1300 Hz during the first ~0.3—1.0 seconds of the 
post-bounce evolution. This is most likely due to an 
avoided crossing [88] of interfering ( = 2 PNS pulsation 
modes that are coupled and mixed [51]. The best cur- 
rent thinking is that the interfering modes are a trapped 


M/M, 
3 We define the compactness here as Êm = Rl [Mo where 


MM) /1000 km’ 
M =1.75 Mo. 


g-mode and the ¢ = 2 f-mode [51, 86], though much work 
remains to be done to determine the details and the na- 
ture of the couplings. Non-linear mode coupling may be 
involved. There is also evidence for some GW power in 
the “bumped” g-mode (that thereafter trends to lower 
frequencies), seen just after the mode repulsion in Fig- 
ures 5 and 6, but most clearly in Figure 17 below. This 
is a qualitatively similar feature to that highlighted re- 
cently in [89]. Nevertheless, within ~0.3 to ~1.0 seconds 
(depending upon the progenitor) most power is clearly in 
the f-mode, where it persists thereafter. 

G-modes, however, are generally at low frequencies be- 
low ~500 Hz and don’t contribute much to the GW sig- 
nature. Moreover, the presence of lepton-gradient-driven 
PNS convection [32, 55] introduces a region in the PNS 
for which g-modes are evanescent and non-propagating. 
Figure 12 depicts the evolution of PNS convection for 
most of the models presented in this work. One sees 
clearly that the extent of PNS convection starts in a nar- 
row shell, but grows wider with time. By ~1.6 to ~4.0 
seconds PNS convection has grown to encompass the cen- 
ter and most of the residual PNS and will persist beyond 
the simulation times of this study [56]. During the early 
post-bounce phase, though most g-modes have frequen- 
cies too low to couple with the f- and p-modes, at an early 
stage before the region of PNS convection has grown too 
thick it is possible for a g-mode trapped mostly interior 
to PNS convection to couple with them. With time, the 
coupling will be broken by the evolving thickness of the 
convective shell; it is this growth that eventually severs 
the coupling with the outer regions where the imping- 
ing plumes are providing the excitation and that leads 
to the jump to the pure f-mode. However, much work 
remains to be done to fully demonstrate the details of 
this coupling and “bumping” transition. Nevertheless, 
the manifest presence of this avoided crossing in the GW 
spectrograms and in the GW signature of core-collapse 
universally is an interesting direct marker of the presence 
of PNS convection that deserves further study. 


E. The Angular Anisotropy of the Matter 
Gravitational-Wave Emissions 


In Figure 13, we plot the matter strain for the 9b, 
12.25-, and 23-Mọ models as a function of time af- 
ter bounce to illustrate the anisotropy with various, 
arbitrarily chosen viewing angles. Note that signifi- 
cant anisotropy manifests at late times in the (low- 
frequency) memory component, which captures the large- 
scale asymmetry in the explosion ejecta. The non- 
exploding 12.25-M, model, by comparison, has virtually 
no anisotropy. The early high-frequency component is 
stochastic, whereas the low-frequency late-time memory 
show secular time-evolution that does not average to zero 
and indicates a metric shift, reaching values of ~5 cm 
for the various massive models (showing a general trend 
with the progenitor mass/explosion asymmetry). We 


find similar significant anisotropy in the (low-frequency) 
neutrino memory component (discussed in §IIF). Im- 
portantly, however, when calculating the total inferred 
“isotropically-equivalent” radiated GW energy, which is 
dominated by the higher-frequency component in and 
near the LIGO band, as a function of angle we find that it 
varies by ~10 to ~15% around an angle-averaged mean. 
This implies that, though the higher-frequency emissions 
are indeed anisotropic, the integrated high-frequency sig- 
nals are only weakly dependent on angle. 


F. The Neutrino Memory Component 


In addition to the matter memory, asymmetries in 
the emission of neutrinos produce another low-frequency 
memory component, the neutrino memory [60-63, 90, 
91]. In Figure 14, we plot the gravitational wave strain 
due to anisotropic neutrino emission as a function of time 
after bounce for the models studied here. The neutrino 
strain is significantly larger in magnitude than the mat- 
ter contribution, reaching over 1000 cm for the most mas- 
sive progenitors. There is generally a hierarchy of strain 
amplitude with progenitor mass, reflecting the sustained 
turbulent accretion in more massive progenitors, which 
results in higher neutrino luminosities and generally more 
anisotropic explosions. The 11-Mọ model is an excep- 
tion and fields the highest strain amplitude. In addition, 
the neutrino memory shows much lower frequency evo- 
lution and more secular time-evolution than the matter 
component, which is fundamentally because it is a cumu- 
lative time-integral of the anisotropy-weighted neutrino 
luminosity (see Appendix A as well as [72]). The differ- 
ence in the mean frequencies of the neutrino and matter 
memories may, therefore, provide a means someday to 
distinguish them observationally. 

In Figure 15, we plot the gravitational wave energy due 
to neutrino emissions as a function of time after bounce 
for our various 3D models. There are several key and dis- 
tinguishable features. First, like the gravitational wave 
energy due to matter motions displayed in Figure 10, we 
see growth by over two orders of magnitude (but not 
three, as in the former) in the neutrino memory. Ad- 
ditionally, we generally see a hierarchy with progenitor 
mass, however with the 11-Mọ surpassing the 23-Mo 
until ~4.5 seconds. Unlike the matter component of the 
gravitational wave energy, the neutrino component shows 
sustained growth for our longest duration model (the 23- 
Mg model). In comparison with Figure 5 from [92], this 
emphasizes again the need to carry simulations out to 
late times to capture the entire signal. Note that, despite 
the higher strains seen in the neutrino component of the 
gravitational wave signature and its sustained growth, 
due to the much smaller frequencies it is still more than 
two orders of magnitude less energetic than the matter- 
sourced gravitational wave energy. In addition, though 
both components capture the development of turbulence, 
the neutrino component does not show a prompt convec- 


tive phase and begins to develop ~100 to 200 ms later 
than the matter component. 

As with the matter component, the neutrino compo- 
nent is most pronounced for delayed explosions of models 
with higher compactness reflecting their more vigorous 
turbulent accretion history and more anisotropic explo- 
sions. 


G. Turbulent Accretion Excites Gravitational 
Wave Emission 


The major excitation mechanism of gravitational 
waves from core-collapse supernovae and black hole for- 
mers is the pounding accretion onto the PNS core 
[41, 43, 48, 49, 51, 52]. As shown in §I1 A, much of the 
GW power comes out at the frequencies associated with 
the pulsational modes of the PNS. To demonstrate the 
correlation of the gravitational strain with the matter ac- 
cretion, we first remove the matter memory by applying 
a high-pass Butterworth filter below 15 Hz to the strains. 
Then, in Figure 16 we plot the turbulent hydrodynamic 
luminosity evolution interior to the shock with this fil- 
tered strain timeline. 

The turbulent hydrodynamic flux is defined following 
[35] as 


1 r 
Bean = ( Gaz +u+ p) Uburo) (2) 


dQ’ 


including the turbulent kinetic energy, the internal en- 
ergy u, and the pressure p. The turbulent velocity is 
defined here as the radial component of the turbulent 
velocity, 


Yonu = ((0" — ue) dae > (3) 
where v” 


rye iS the density-weighted angle average of the 
radial velocity. This is calculated at 110 km for all models 
except the non-exploding 12.25- and 14-Mo progenitors, 
whose shocks early on sink below this radius; for these 
models we calculate the turbulent hydrodynamic flux at 
2.5 times the PNS radius (defined here as the density 
cutoff at 101° g cm~’). 

The strong correlation throughout their evolution 
(even in detail) between the turbulent hydrodynamic flux 
impinging onto the core and the GW strains demon- 
strates that turbulent accretion through the shock and 
onto the PNS core is the major agency of GW excitation 
and emission in core-collapse supernovae. We note that 
this correlation was demonstrated even though the flux 
was angle-averaged and the strain was for emission along 
the x-axis. No attempt was made to break the flux into 
L = 2;m = |—2, —1,0, 1,2] components, yet the correla- 
tion is clear. 

To provide another perspective, in Figure 17, we show 
the relation between the gravitational wave energy spec- 
trogram and the accretion rate “power” in frequency 


components with f > 25 Hz (orange lines) for the 9.5- 
and 23-Mọ models. As the panels demonstrate, after the 
explosion and the period of heavy infall subsides (dur- 
ing which the effects of individual accretion events over- 
lap), the accretion rate power shows clear correlations 
with excursions in the gravitational wave energy spec- 
trogram. Spikes and gaps on the spectrogram coincide 
directly with the peaks and troughs on the curve, mean- 
ing that the gravitational waves are excited by, or at least 
correlated with, the short-period variations in accretion 
rate onto the core. Such mass accretion rate variations 
directly tie both episodic fallback and outflow events with 
the GW emission. 

We note that the colormap used for these plots re- 
veals, particularly for the 23-Mọ model, some power in 
the g-mode “bumped” by the f-mode and the repulsion 
between the two modes around ~0.8 seconds (see §IID). 
Though weak, for the 23-Mo model this signature con- 
tinues almost to ~2.0 seconds, and perhaps beyond. 


1. Possible Secondary Role of Proto-Neutron-Star 
Convection 


Some have suggested that inner PNS convection itself 
excites much of the GW emission from CCSNe [47]. In 
Figures 18 and 19, we overplot the angle-averaged con- 
vective hydrodynamic luminosity at its peak inside the 
PNS (see also Figures 12) against the Butterworth fil- 
tered strain. We follow Eq.3, but limit the radius to 
that of the maximum convective luminosity within the 
PNS (positive outwards), which lies near ~15 km for the 
various models. As Figures 18 and 19 show there seems 
to be no correlation between the two. This is particu- 
larly clear at late times when the PNS convective flux 
is still large while the GW emissions have all but sub- 
sided. If PNS convection were the agency of excitation 
at all phases, GW emission would not have subsided to 
such a degree at later times (from 0.3 to 1.5 seconds after 
bounce). This comparison demonstrates the importance 
of simulating to late times to capture the entire GW sig- 
nal. 

However, as Figures 5 and 6 themselves show, the f- 
mode persists to the latest times and manifests an episod- 
ically modulated (see Figure 17), though continuous, sig- 
nal. We have yet to identify the major excitation mecha- 
nism for this component at late times. Anisotropic winds 
and neutrino emissions from the residual core could be 
causes, but inner PNS convection may be also a factor 
here. It is also the case that the mode should ring down 
over a period given by the dominant damping mechanism. 
Such mechanisms include sound generation, the back re- 
action of anisotropic winds, neutrino emission coupling 
and viscosity, non-linear parent-daughter mode coupling, 
and numerical dissipation. We reiterate, however, that 
the f-mode continues to ring and produce a weak GW 
signal for the duration of all our simulations. Clearly, 
this topic deserves more detailed scrutiny in the future. 


Nevertheless, this later phase amounts to only a few per- 
cent of the total energy emitted. 


III. CONCLUSIONS 


In this paper, we have presented and analyzed the 
gravitational-wave signatures of an extensive suite of de- 
tailed initially non-rotating 3D core-collapse supernova 
simulations spanning a wide range of massive-star pro- 
genitors. For the first time, most of the published sim- 
ulations were carried out to late enough times to cap- 
ture more than 99% of the total GW emission from 
such events. Moreover, we have endeavored to dump 
the relevant quadrupole data at a rate sufficient to ef- 
fectively eliminate Nyquist sampling errors. We see that 
the f/g-mode and f-mode oscillation modes of the PNS 
core carry away most of the GW power and that generi- 
cally there are avoided crossings and modal interactions 
likely associated with the evolution, extent, and charac- 
ter of lepton-driven PNS convection. The f-mode fre- 
quency inexorably rises as the proto-neutron star core 
shrinks during its Kelvin-Helmholtz contraction phase, 
driven by neutrino loses, and its power and frequency 
behavior are central features of the GW emissions from 
the core-collapse event. Other modes are also seen in 
the GW spectra, in particular a £ = 2;n = 1 p-mode 
and, perhaps directly, a trapped g-mode, though most 
g-modes are not excited. Whether other p-modes are in 
evidence is to be determined. 

We demonstrate that the GW emission is powered 
mostly by accretion plumes onto the PNS that excite its 
modal oscillations and also produce a “haze” of higher 
frequency emission also correlated with the phase of vio- 
lent accretion, after which the signal subsides to be domi- 
nated by the chirp of the f-mode signal at low power that 
nevertheless continues beyond the duration of even these 
simulations, albeit weakly. The duration of the major 
phase of emission varies with exploding progenitor and 
is generally shorter for the lower-mass progenitors (~0.3- 
0.5 seconds) and longer for the higher-mass progenitors 
(~1.5 seconds). We find a strong correlation between the 
total GW energy radiated and the compactness of the 
progenitor whose mantle explodes as a supernova. Fur- 
thermore, we find that the total GW energy emissions can 
vary by as much as three orders of magnitude from star 
to star. Hence, there is a severe progenitor dependence 
that must be factored into any discussion of detectabil- 
ity. For the black-hole forming models, since accretion is 
not reversed at any time or at any solid angle, their GW 
signal lasts until the black hole forms, tapering off only 
slowly until then. In addition, they do not manifest the 
high-frequency haze seen for the exploding models. For 
these black-hole formers, we also witness the emergence 
of a spiral shock motion that modulates the GW emis- 
sion at a frequency near ~100 Hz that slowly increases 
as the stalled shock sinks. 

In Figure 20, we plot the sensitivity curves with the 


amplitude spectral densities at 10 kiloparsecs of our 3D 
models. More massive models generally leave larger 
footprints, with the 11 Mo studied being the excep- 
tion. Current and next-generation detectors can ob- 
serve, for galactic events, a signature spanning orders of 
magnitude, from subHz to ~3000 Hz. While Advanced 
LIGO/Virgo/Kagra (LVK, [93]) can detect ~30—3000 
Hz signals for the more massive progenitors, upcoming 
detectors, including the Einstein Telescope [94, 95], the 
Cosmic Explorer [96], BBO [97], and Decigo [98, 99] (at 
lower frequencies, see also [100]) should be able to de- 
tect galactic events for all progenitor masses studied here 
through almost three orders of magnitude in both fre- 
quency and total energy radiated. However, a detailed 
retrieval analysis, informed by the best signal-processing 
approaches, has yet to be performed and is an important 
topic for future work. 

Though we have endeavored here to provide a compre- 
hensive look at the gravitational-wave signatures of core 
collapse, there remains much yet to understand. Topics 
unaddressed here are the nuclear equation-of-state de- 
pendencies, the role of rapid rotation, the possible sig- 
natures of strong magnetic fields [36, 89](REF), and the 
results for other progenitor massive stars and from other 
stellar evolution codes. Importantly, the analysis of a de- 
tected GW signal would be significantly aided if done in 
concert with a corresponding analysis of the simultaneous 
neutrino signal. Optimal methodologies with which to 
extract physical information from such an analysis have 
yet to be designed. Nevertheless, we have in the GW 
signal of core-collapse supernovae a direct and real-time 
window into the supernova mechanism and PNS evolu- 
tion. Therefore, such a methodology would likely pay rich 
scientific dividends when astronomy is finally presented 
with the opportunity to employ it. 


DATA AVAILABILITY 


The numerical data associated with this article and 
sonifications of the gravitational wave strains will be 
shared upon reasonable request to the corresponding 
author. The gravitational wave strains as well as 
the quadrupole data are available publicly at https: 
//dvartany.github.io/data/. 
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Appendix A: General Equations 


We follow the method in ([76], see also [47, 72, 77, 102]) 
to calculate the gravitational wave strain tensor. We cal- 
culate the first time derivative of the mass quadrupole 
with the following formula: 


F 2 
a) = ŽQy= fè (væg tozi Sor), (AL) 


where v is the velocity, x is the Cartesian coordinate, 
Qij is the transverse-traceless quadrupole tensor, and r 
is the radius. The transverse-traceless gravitational wave 
strain tensor hg? is calculated by taking the numerical 
time derivative of qij, i.e., 

HEP ws 2G dqij 


IT — A2 
ij éD dt’ ( ) 


where D is the distance to the source. Hereafter, we 
drop the superscript “TT”. We also calculate and dump 
the quadrupole Qi; = f dx p(xixj— år?6;j), and its nu- 
merical derivatives are consistent with the values calcu- 
lated by the above equation. However, taking numerical 
derivatives can be viewed as convolving the signal with 


a window function, thus it will introduce a bit of extra 
noise at high frequencies. 

The “plus” and “cross” polarized strains along (0, 6) 
direction are given by 


G dqee ddog 
hy = A 
ii a5 ( dt dt oe) 
and 
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where (from [77]) 


qoo = (qzx CoS? $ + dyy sin? 6 + 2 dey Sin ¢ cos 6) cos? 0 


+ qz- sin? 0 — 2 (qxz CoS d + qyz Sin œ) sin 6 cos 0 
(A5) 


dop = qax sin? 6 + dyy cos? p— 2 dxysinpcosøp (A6) 


qog = (dyy — Gea) cos O sin ġ cos ø + dry cos O (cos? ¢ — sin? ø) 


+ dxz Sin 0 sin @ — qyz sind cos @. 
(A7) 


The total energy emitted in gravitational waves is 
t 3 2 
G (Pij 
Ecw = 7 
m I 2 Bes ( dt? ) (A8) 


To calculate gravitational waves from neutrino asym- 
metries, we follow the prescription of ([91], see also [92]). 
We include angle-dependence of the observer through the 
viewing angles a € [—7,7] and 8 € [0,7]. The time- 
dependent neutrino emission anisotropy parameter for 
each polarization is defined as 


aslt,a, 8) = fy fa! Ws(O',0,8) (Mt), 


(A9) 
where the subscript is S € {+, x} and the gravitational 
wave strain from neutrinos is defined as: 


t 
hs(t,0,8) = FÉ I dt’ A(t) as(t,a, B), (A10) 


where A(t) is the angle-integrated neutrino luminosity as 
a function of time, D is the distance to the source, and 
Ws (’,a, 3) is the geometric weight for the anisotropy 
parameter given by 


Dg(6", ¢', a, B) (A11) 
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where (from [91]) 


D, = [1 + (cos(¢’) cos(a) + sin(¢’) sin(a@)) sin(@’) sin(8) + cos(6’) 
cos(8)]{[(cos(¢’) cos(a) + sin(¢’) sin(a)) sin(’) cos(8) 
— cos(6’) sin(8)]? — sin?(’)(sin(¢’) cos(a) — cos(¢") sin(a))?} 


(Al2a) 


D, = [1 + (cos(¢’) cos(a) + sin(¢’) sin(a)) sin(@’) sin(8) + cos(6’) 
cos(8)}2[(cos(¢’) cos(a) + sin(¢’) sin(a)) sin 
— cos(6’) sin(8)] sin(9’)(sin(¢’) cos(a) — cos 


6") cos(8) 
VW) sin(a))” 


= [(cos(¢’) cos(a) + sin(¢’) sin(a)) sin(6’) cos(8) — cos(0") 
sin(3)]? + sin? (6’)(sin(¢’) cos(a) — cos(¢’) sin(a))? . 


(A12c) 
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TABLE I. Model Properties 
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Model Run Time Explosion? Shock Velocity  f/g-mode Nyquist 
Energy Fraction Frequency 
(Mo) (s, pb) (km s~') Total (Late) (Hz) 
90 1.47 v 16,000 56.08% (N/A) 5000 
9BW 1.1 v 14,000 68.11% (N/A) 8000 
9F 1.6 v 16,000 60.92% (5.75%) 3000 
9F,1 PAR v 16,000 76.68% (6.63%) 5000 
9.25 2.7 v 13,000 72.55% (3.06%) 8000 
9.5 2.1 v 11,000 76.35% (0.72%) 8000 
11 4.5 v 11,000 85.02% (13.55%) 6000 
12.25 2.0 x 94.55% (25.71%) 8000 
14 1.0 x 75.82% (N/A) 8000 
15.01 2.0 v 6,000 84.46% (5.12%) 8000 
23 6.2 v 6,000 85.74% (8.50%) 8000 


Table of some the features of the 3D simulation. We include the simulation time, in seconds 
post-bounce, the explosion outcome, the asymptotic shock velocity, and the fraction of the GW 
energy radiated via the f/g-mode. Models with a checkmark explode, and models with an X do not 
explode. The various 9 solar mass models were run on Theta (ALCC), Blue Waters (NCSA), and 
Frontera (TACC), respectively. 9F,] is a longer simulation of the 9-Mo progenitor, also done on 
Frontera. The f/g-mode energy is the gravitational wave energy within +100 Hz of its central 
frequency. “Total” means for all time, while “Late” means after 1.5 seconds after bounce. After 
the mode repulsion, almost all the power is the an unmixed f-mode. Also in the table are the 
lowest Nyquist sampling frequencies for each run. The behavior at frequencies below these should 
not be compromised by Nyquist errors. 
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FIG. 1. Mass density profiles for the various progenitors studied here, from 9.0 (reddish) to 23 (blue) Me. The accretion of the 
Si/O interface frequently corresponds with the onset of shock revival. Note that the 12.25 and 14 Me models do not explode 
and are destined to produce a stellar-mass black hole. 
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FIG. 2. Early post-bounce behavior of the mean shock radii (in km) for all the models as a function of time (in seconds after 
core bounce). All models explode except for the 12.25- and 14-Mọ progenitors, which show evidence for a spiral SASI after 
~400 and ~300 ms post-bounce, respectively, in the form of ~10 ms oscillations in the shock radii. Note that the 12.25- and- 


Mo (non-exploding) progenitors also shows a secular ~70 ms oscillation in the shock radii. 


17 


h, x x D [cm] 


0 1 2 


3 4 
Time after bounce [s] 


FIG. 3. The gravitational wave strain, from matter motions, times distance (product in cm) for all the 3D models studied in 
this paper as a function of time after bounce (in seconds) observed along the positive x-axis in the simulation frame. Both cross 
(blue) and plus (black) polarizations are illustrated. The red star indicates the approximate time of the onset of explosion, 
defined here as when the stalled shock radius changes concavity at ~150 km and accelerates. See Figure 2. 
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FIG. 4. Same as Figure3, but now focusing on the matter strain during the first 250 ms after core bounce to illustrate the 
prompt convection. Note the stately oscillation present in most models of ~10 ms periods visible in the + polarization and 
absent in the x polarization. This is a generic feature indicative of the strain polarization geometry. Note that only model 9a 
has any perturbations on bounce. 
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FIG. 5. In this figure and the next, we provide the gravitational wave energy spectrogram (to a normalization constant of 
order unity) at 10 kiloparsecs (kpc) as a function of time after bounce (in seconds) for the 3D models highlighted in this 
paper. Note that the frequency ranges and durations differ from panel to panel. The colormap provides the spectral density 
in units of Bethes (10°! egs) per Hertz, where we have integrated over a running 40-millisecond bin around each time. We 
have found that the results do not depend significantly upon the width of this window. Though the models in this first set 
are for 9-, 9.25-, and 9.5 Mo, the description in this caption also applies to the 11-, 12.25-, 15.01-, and 23-Mọ models in the 
accompanying figure (Figure 6). All models show a strong feature out to ~50 ms and up to ~2000 Hz associated with prompt 
convection, followed by a short quiescent phase, during which the turbulence between the shock and the proto-neutron star 
grows. Thereafter, all models illustrate high frequency behavior before ~one second (even earlier for the rapidly exploding 
lower progenitor mass models), which truncates in power as the turbulent accretion which excites it subsides. This “emission 
haze” is either a superposition of core pulsational £ = 2 p-modes with a range of node numbers or the GW emission from 
the supersonic plumes themselves impinging upon the PNS core, or some combination of both. The fundamental f/g-mode 
(early), then the f-mode (later), is manifest and significant in all models for the duration of the simulations. We note that 
the 15.01-Mo model (next figure), which has a relatively lower ejecta mass, sustains turbulent accretion and high frequency 
haze emission until the end of our simulation, ~2.0 seconds post-bounce. We also note that we see higher harmonics/overtones 
after the first second and after the earlier emission haze abates, most visibly in the 15.01- and 23-Mo models. These later-time 
PNS pulsational modes are likely the £ = 2,n = 1 p-mode. The fundamental f-mode frequency grows roughly quadratically 
with time ([102]) and is associated with the shrinkage of the PNS as it executes its Kelvin-Helmholtz contraction phase, driven 
by neutrino losses. Except for the very early phase after bounce, most g-modes are suppressed due both to the presence and 
growth of PNS convection in the deep core and to their lower frequencies (~a few hundred Hz). Early on, all models show a 
dark band near ~1000 Hz which may be a signature of an avoided crossing and interference with a trapped g-mode. See text 
in §II A and §IID for a discussion. 
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FIG. 6. Same as Figure 5, but for models 11-, 12.25-, 15.01-, and 23-Mo. 
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FIG. 7. Plots (2D slices of 3D data) of the Mach number distribution in the inner 100 kilometers of the 23-Mọ simulation at 
574 and 632 milliseconds after bounce. Seen are representative supersonic plumes crashing onto the PNS core. Such plumes 
not only excite the GW emission, but themselves emit a component of the GW signal. See text in $II A for a discussion. 
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FIG. 8. Plots of the effective strains as a function of time after bounce (in seconds) for two representative 3D models (the 
black-hole former 14-Me and late-exploding 23-Mo model). Note that we see a dark red band at low (<50 Hz) frequencies 
for the exploding model associated with the asymmetric ejecta at later times (matter “memory”). The non-exploding model 
shows a spiral SASI signature at ~100 Hz. This quasi-periodic feature is also seen after ~0.4 seconds in Figure 2 for both the 
14-Mo and 12.25 Mọ black-hole formers. Such a feature shows up only for non-exploding models when the shock radius has 
shrunken below ~100 km. These features are much less prominent in the gravitational wave energy spectra because they are 
low-frequency, and, hence, low-power components. There may be a hint of the SASI itself in the lower left-hand corner of the 
panel to the right. See text in §II A for a discussion. 
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FIG. 9. In this figure, we show the spectrogram of the z-component of shock dipole position for four models. Note the focus 
on lower frequencies. In the non-exploding 12.25 and 14 Meo models, the spiral-SASI modes are clearly seen and manifest 
gradually increasing frequencies as the PNS core shrinks. The exploding models (for example, the 15.01 and 23 Me models 
depicted here) don’t show this feature, but do manifest the even lower frequency signals due to the matter memory caused by 
explosion anisotropy (see also Figure 8). We note that the 23-Mo model may show a weak SASI signature (analogous with that 
found in the 25-Mo model in reference [72]) that vanishes at ~ 0.55 seconds and is perhaps indicative of its delayed explosion 
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FIG. 10. The gravitational wave energy from matter motions (in Mọ c°) for all the 3D models highlighted in this paper as 
a function of time after bounce (in seconds). Note that the total gravitational wave energy radiated differs by ~three orders 
of magnitude from the least massive, 9-Mo to the most massive, 23-Mọ progenitor. This energy grows by three orders of 
magnitude for the most massive progenitors over the first ~2 seconds of simulation, but is already asymptoting shortly after 
one second post-bounce. All models show rapid growth in the first ~50 ms, associated with the onset of prompt convection 
driven by the overturn of the shocked mantle dynamically generated at and after bounce as the shock stalls initially into 
accretion. After this phase, the neutrino-driven turbulence between the shock and the proto-neutron star core grows in vigor 
over a period of ~100 milliseconds and excites a spectrum of core pulsational f- and some p-modes and likely generates a GW 
component due to the impinging of the plumes onto the PNS that all together constitute the bulk of the gravitational radiation 
issuing from the supernova. This phase can last from hundreds of milliseconds to <1.5 seconds, depending upon progenitor, 
after which the strains subside to a hum dominated by the fundamental £ = 2 f-mode and (weakly) overtones. This last phase 
can last for many seconds. However, the signals from those models destined to leave black holes are still vigorous for a longer 
period of time, since accretion persists for these models until the black hole forms, after which the signal ceases abruptly. See 
the text for a discussion. 
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FIG. 11. The gravitational wave energy from matter motions (in Moc?) for all the 3D models in this study versus the model 
compactness at 1.75 Me [74]. Note the clear linear correlation on the log-log scale, with a power-law slope of ~0.73. The outliers 
are the non-exploding models (12.25- and 14-Mo, shown with triangle markers) and these models generally take longer to radiate 
a given amount of gravitational-wave energy (see Figure 10). The black-hole formers don’t as a rule manifest turbulence that 
is as vigorous as those models that eventually explode. All exploding models show a steep growth in the gravitational-wave 
energy associated with the onset of vigorous turbulence. The 23-Mo, which explodes the latest, has the latest growth spurt 
in gravitational energy. By contrast, non-exploding models show a slower rate of of increase of gravitational-wave emission, 
associated with a slower growth rate of turbulence. However, due to the persistently stalled shock and continued accretion, the 
black-hole formers will generate gravitational waves for a longer time and will saturate their total gravitational-wave emission 
more slowly. Finally, a spiral SASI mode emerges for the black-hole formers when the stalled shock wave retreats significantly 
(below ~100 km). This mode modulates the stalled accretion shock radius and the gravitational-wave emissions of these 
black-hole progenitors quasi-periodically. 


26 


10! 


10° i me — 10° 
107 — 107 = 
n $ y = n 
107? 30 = 107 20 
oO g ú 2 
—10°2 2 —10°? $ 
3 3 
joe - in —107 4 
—10° . | —10° 
10! =, i 7 10! 
10 i : T : oa 10 
t — thounce [S] t — thounce [S] 
89.5 i s11.0 Å 
10 TL nL EE A GE EES RE FEEL EC 10 
10° te 10° 
lot = ; " a ‘ 107 = 
= 102? mp = 1.25 F ø w? 4 
T 0 8 ~ 1.00 H 0 3 
8 [=] [=] 
= -10 > = Oo p -10 2 
3 3 
-10 4 50 Hate gt eN 
10° .25 ARI —10° 
10! i i i i 40! 
10 ) 5 i F - 10 
t — thounce [S] 
s14.0 
10! 2.00 peg oe 10! 
10° 7 10° 
0" = . wi 
> 2 > 2 
10? 2p = 1.2! ? 10? spg 
0 a A " 0 a 
© = 
a a 
—10” + + 0. niii a i —10” > 
-107 4 J ii baht 107 = 
—10° ,2: i Wa —10° 
—10! 0 0.00 Eee —10! 
t — thounce [S] t — thounce [S] 
815.01 ; s23.0 i 
2.00 mmr 10 15 et 10 
1.75 10° 10° 
1.50 io 5 
» 2 
10° so 
l] 
0 g 
oO 
a 
-10> + 
3 
0.50 —1907 2 
0.25 i —10? 
0.00 Pi Pi ta (aaa Eai —10! 1 AL aif LA 0 4 —10! 


it 
0.50 0.75 1.00 1.25 1.50 i 2 3 
t — thounce [8] t — thounce [S] 


FIG. 12. Spacetime diagrams illustrating the convective regions in the PNS for eight representative models, identified by the 
convective luminosity, as a function of time after bounce. Note that convection reaches the PNS core within ~1.6 to ~4.0 
seconds, earlier for the low-mass progenitors and later for the 23-Mọ model. For all models, PNS convection starts in a shell 
and then grows to encompass much of the core. The fact that g-modes are evanescent in convective regions bears on the 
suppression of g-mode excitation in general, but also on the eventual severing of the interaction/coupling between a trapped 
g-mode interior to it and the PNS surface where the plumes are exciting much of the GW emission. See text for a discussion 
(SII D): 
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FIG. 13. The gravitational wave strain h;D (cm) as a function of time after bounce (in seconds) for the 9b model, the non- 
exploding model 12.25-Mo, and model 23-Mo@ for various, arbitrary, viewing angles illustrated by color. These plots do not 
include the neutrino memory (§II F). Note the large anisotropy by viewing angle, which is dominated at late times by the matter 
memory, which is sourced by the large-scale explosion asymmetry. In the black hole model (12.25-Mo), we see no memory, and 
little anisotropy, at late times. Though we show only the + polarization, the results are consistent for both polarizations. We 
see that in principle some measure of the degree and angle of the matter explosion asymmetry can be discerned by the sign and 
magnitude of this memory. However, the neutrino memory, if only weakly correlated with the matter memory, will complicate 
the extraction of such a feature. Nevertheless, the difference in the frequency spectra of the matter and neutrino memories 
generically, with the former at a higher average frequency, may allow one someday to disentangle the two. 


28 


1200- T T | 
1000 
800 
600 


400 


+x D [cm] 


= L l l l l l l 
4005 1 2 3 4 5 6 


Time after bounce [s] 


QL 


FIG. 14. Same as Figure3, but for the gravitational wave memory due to neutrino emission anisotropy. Because of the 
cumulative nature of the neutrino gravitational wave memory, the strain is both larger and shows less variation on small 
timescales when compared to the matter-driven strain. We reach strains as high as 1000 cm for the 11- and 23-Mọ models. 
Note that this is consistent with earlier work [92] showing continued growth with time of the memory strain. The strain growth 
rate also accelerates, showing upwards of 600 cm growth from ~2 to ~3 seconds for the 11-Mo progenitor. This emphasizes 
our essential point that models should be carried out to late times to fully capture not only the explosion dynamics, but the 
emerging gravitational wave signatures. The 11-Mo progenitor has the highest neutrino memory strain among those models 
we present in this paper. Otherwise, we see a general correlation with mass, with the low-mass progenitors having the smallest 


neutrino memory. The non-exploding 12.25- and 14-Mọ models show significant evolution and magnitude in their memory 
signatures. 
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FIG. 15. Same as Figure 10, but for the radiated gravitational wave energy due to the neutrino anisotropy. As seen in Figure 14, 
unlike the matter contribution, the neutrino anisotropy continues to contribute to the gravitational signal past six seconds (for 
our longest simulation), emphasizing the presence of sustained neutrino emission asymmetry. We see again at least two orders 
of magnitude of growth with time after —50 ms in the associated gravitational wave energy, this time in the neutrino memory 
component. The 23-Mo progenitor has neutrino GW memory energies only slightly higher than the matter GW components 
for the least energetic, lowest-mass progenitors. Importantly, for the same progenitor model, the late time neutrino-sourced 
gravitational-wave energies, albeit still growing, are more than two orders of magnitude smaller than the corresponding matter 
components. 
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FIG. 16. Turbulent luminosity (blue, in 10° erg s7") at 110 km (except for model 12.25-Mo, where it is evaluated at 2.5 times 
the PNS radius), overplotted with both filtered gravitational wave strain polarizations (green and black) multiplied by distance 
(product in cm) as a function of time after bounce (in seconds). We emphasize the clear correlation shown here — the outer 
turbulent luminosity closely follows the strain. It reproduces both the global time behavior, rising and following at coincident 
times with the strain, and the episodic ‘packets’ of accretion (for instance, reproducing the bump in the 15.01-Mo strain at 
~1.35 s and for the 23-Mo at ~3.5 s). From this figure, it is clear that it is the outer downward turbulent convective flux onto 
the PNS that drives the gravitational wave signature after the prompt convection phase. 
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FIG. 17. In this figure, we show the correlation between the gravitational wave energy spectrogram and the accretion rate 
power in frequency components with f > 25 Hz for two representative models. The background spectra in this figure are the 
same as those plotted in Figure 5 and 6, but with different color-bar ranges so that the late time structures can be more clearly 
seen. The orange curve shows the accretion rate “power” in frequency components with f > 25 Hz. The 25-Hz filter is chosen 
here because the spectrogram is calculated using a 0.04s-wide sliding window. The accretion rates are measured at a radius of 
25 km, and the PNS radii of the two models fall below this value at ~0.9 seconds and ~1.0 seconds after bounce, respectively. 
Note also that this colormap for the 23-Mo model spectrogram reveals a weak signal below the f-mode associated with the 
avoided crossing that may be the bumped g-mode and that seems to continue beyond ~1.5 seconds. See text in §IIG for a 
discussion. 
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FIG. 18. Convective hydrodynamic luminosity within the PNS (blue, in 10°! erg s~') at its peak value in the inner PNS, 
overplotted with both gravitational wave strain polarizations (black and green) multiplied by distance (product in cm) as a 
function of time after bounce (in seconds) for models 9a, 9b, 9.25, and 9.5. Figure 19 depicts the same quantities, but for models 
11-, 12.25-, 15.01-, and 23-Mọ. Note that PNS convection first increases in strength over time and then stays approximately 
constant (with some minor excursions) or stays high with a secular evolution for the duration of the simulations. Moreover, 
PNS convection grows in extent over time to encompass more and more of the core, expanding deeper in radius and reaching 
the center of the PNS after from ~1.7 seconds (the 9-Mọ model) to ~4 seconds (the 23-Mg model, see Figure 19), all the 
while maintaining its vigor. This behavior shows little or no correlation with the gravitational-wave strain, which decreases in 
magnitude with time from an earlier peak and nearly flattens at late time (dominated then by a humming fundamental {mode 
(see Figure 5). In particular, at late times, when the strain has subsided substantially, the inner PNS convective luminosity 
is always still strong, indicating that PNS convection is not the dominant excitation mechanism of CCSN gravitational-wave 
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FIG. 19. Same as Figure 18, but for models 11-, 12.25-, 15.01-, and 23-Mo. 
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FIG. 20. Gravitational wave sensitivity curve for all the models as a function of time (in seconds after core bounce), compared 
with the CCSN GW signals at 10 kpc, including both matter (high frequency, dominant above several hundred Hz) and 
neutrino contributions to gravitational wave emissions (low frequency, dominant from sub-Hz to several tens of Hz). Note that 
DECIGO, the Einstein Telescope, and the Cosmic Explorer can probe down to the lowest-mass progenitors across three decades 
in frequency (~10 to ~5000 Hz), whereas aLIGO is sensitive between ~30 to ~3000 Hz only to the most massive progenitors. 


